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ABSTRACT 

We reinvestigate the problem of the appearance of relativistic jets when geometrical 
opening is taken into account. We propose a new criterion to define apparent veloc- 
ities and Doppler factors, which we think being determined by the brightest zone of 
the jet. We numerically compute the apparent velocity and the Doppler factor of a 
non homokinetic jet using different velocity profiles. We argue that if the motion is 
relativistic, the high superluminal velocities P app ~ 7, expected in the case of an 
homokinetic jet, arc only possible for geometrical collimation smaller than the rela- 
tivistic beaming angle 7 . This is relatively independent of the jet velocity profile. 
For jet collimation angles larger than 7 _1 , the apparent image of the jet will always be 
dominated by parts of the jet traveling directly towards the observer at lorentz factors 
< 7 resulting in maximal apparent velocities smaller than 7. Furthermore, getting 
rid of the homokinetic hypothesis yields a complex relation between the observing 
angle and the Doppler factor, resulting in important consequences for the numerical 
computation of AGN population and unification scheme model. 

Key words: galaxies: active - BL Lacertae objects: individual: - galaxies: jets - 
gamma-rays: theory - radiation mechanisms: non-thermal 



1 INTRODUCTION 

Jets are present in a very wide number of astrophysical 
objects, from the young stellar objects and X-ray binaries, 
at the galactic scale, to the powerful Active Galactic Nuclei 
(hereafter AGN) whose jets can cross hundreds of kpc in the 
intergalactic medium. If the jets are in relativistic motions, 
Doppler and time delay effects will greatly affect their 
appearance. From a taxonomic point of view, the presence 
or absence of jets, and the inclination of the jet with the line 
of sight, is commonly used to classified extragalactic objects 
in different categories (e.g. lUrrv fc Padovanil Il995l : llTrrvl 
l2005t ) even if it is yet not clear that the jet characterizes a 
particular type of AGN or only a period of the cosmological 
evolution of galaxies. 

Relativ istic motion has been proposed very early by 
iReesI i|l966h to solve the mystery of high luminosities and 
very rapid variability of extragalactic radio sources, which 
should have resulted in a catastrophic Inverse Compton 
cooling if the sources were at rest with respect to the 
observer. He predicted then the possibility of observing 
superluminal motions, a prediction beautifully confirmed 
thanks to the development of high resolution Very Large 
Base Interferometry. Later, detection of high energy 
gamma-ray radiation by EGRET led to the conclusion 



that relativistic motion was also necessary to escape the 
problem of self-absorption by pair production. The most 
simple model that reproduces relatively well the observed 
Spectral Energy Distribution (hereafter SED) is the so- 
called one-zone Synchrotron self-Compton (SSC) model. It 
assumes a spherical zone in which high energy particles are 
injected and cooled via synchrotron and Inverse Compton 
processes. The model parameters are the zone radius, the 
magnetic field as well as the density and the distribution 
of the particles. To reproduce the observed SED however, 
a static model is generally ruled out due to causality 
constraints (see lZensuslll997i and references therein) and a 
necessarily low pair creation optical depth is needed in order 
to avoid all gamma-rays p hotons being absorbed to form 
electr o n-positron pairs (e.g. Maraschi. Ghisellini. fc Celottil 



1 19921 : iHenri. Pelletier. fc Rolandl I1993I V The solution is 
then to assume that the source is uniformly moving with 
a relativistic bulk velocity v = /3c. In this case, all specific 
intensities, in the source rest frame, are a factor 8 Z lower 
than in the observer frame, S — I/T(I — Pfi) being the 
Doppler beaming factor, F the usual Lorentz factor and 
/i = cos 6 bs the cosine angle of the jet axis with respect 
to the observer line of sight (LOS). So the actual photon 
density in the jet frame is much lower than what would 
be deduced for a static source. The possible geometrical 
opening of the jet and non-uniformity of the Lorentz factor 
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across the jet section are neglected in the simplest form of 
this model. 

However jet opening angles are observed in different 
types of ob j ects li ke AGN or Young Stellar Objects 
ifjunor et all Il999l ; iHoriuchi et~ai1 \20o4 ). Some of these 
observations show a decrease of the jet opening angle 
with distance from the central core, an indication of some 
collimation processes. Jet models also predict a variation of 
the j e t opening angle ( e.g. lFerreiralll997 : ICasse fc Keppens 



X 



20021 ; lMcKinnevll2006l ; iHawlev fc Krolikl 120061 ; IZanni et a" 



2007h and in deed some of them ar e able to reproduce the 



observations (|Dougados et al.ll2004 ) 

Given the strong dependance of the Doppler beaming 
factor on the angle between the line of sight and the direction 
of displacement of the emitting region, taking into account 
the jet opening angle should have important impacts on the 
obs erved j et emission. In a series of papers, Gopal-Krishna et 
al. (|2007al and reference therein. Hereafter GK) have inves- 
tigated some of these effects on the observed parameters of 
blazar jets like the jet orientation angle, its apparent speed 
and Doppler factor. 

These authors proposed to compute the effective 
apparent velocity of the jet as a Doppler-boosted inten- 
sity weighted average of the appar ent velocity of each 
point of the jet blob (cf Eq. (5) of iGopal- Krishna et al.l 
l2007al ). We think however that this prescription does not 
really reproduce the way by which apparent velocities are 
measured, i.e. by the displacement of the maximum of 
the image fit. Moreover, due to the light travel effects, 
neglected in the previous analyses, the shape of the emit- 
ting region is different between the jet and the observer's 
frame, hence modifying the apparent shape on the sky plane. 

In this paper we propose to re-analyse these different 
effects using our new prescription to estimate the apparent 
jet velocity. For sake of simplicity, we use a simple formal- 
ism for the jet geometry described in Sect. 12.11 and assume 
an infinitely thin shell propagating in the jet, which would 
capture the essential features of any perturbation traveling 
along the jet. We compare four different jet velocity profiles 
to see the influence of precise angular distribution of Lorentz 
factors. We detail the way we compute the apparent velocity 
of the whole pattern seen by an observer at infinity in Sect. 
12.21 We then present our results in Sect. [3] and discuss them 
in Sect. [4] before giving concluding remarks in Sect. [5] 



2 THE MODEL 
2.1 Formalism 

We consider, in the following work, the simple case of a 
instantaneous central perturbation, propagating with a 
relativistic speed characterized by the Lorentz factor 70 
on the jet axis. In the jet rest frame, we assume that the 
surface emissivity is uniform, its spectrum being described 
by a power-law of index n defined by S v oc v~ n (see 
Sect. 1531 for the case of an angular dependent emissivity). 
The geometrical collimation of the jet is characterized 
by a characteristic parameter 9j, whose exact definition 




Observer 



Figure 1. Sketch of the jet model in the case of velocity distribu- 
tion D\. See text for the signification of the different parameters. 



depends on the assumed shape. The velocity distribution 
is described by a function 7(0), where is the angle to the 
jet axis. For a point of this surface at a given polar angle 0, 
the velocity vector also points in the radial direction (see 
Fig. E]). 

Four different velocity profiles will be used to study 
the influence of the jet structure on the observational pa- 
rameters, and compared with the single blob, homokinetic 
results. The first one (-Di) considers the simple case of a 
conical surface with a constant Lorentz factor 70 inside a 
cone of half-opening angle 9j , and a Lorentz factor equal to 
unity outside (i.e. null velocity). This distribution is obvi- 
ously not physical due to the sharp discontinuity at the edge 
of the cone, but is convenient to use. It is the same as the 
one used in most gamma-ray bursts models. In the absence 
of any specific form argued for in the literature, we have also 
studied different, smoother distributions. The second distri- 
bution (D2) is built using the conical distribution Di, and 
adding a power-law decrease of index —2 for the Lorentz fac- 
tor outside the cone. The two last velocity distributions are 
gaussian (-D3) and lorentzian (-D4) profiles. Then for every 
6*| < 7r/2 (approaching jet), we define: 



Di : 7(9) 
D 2 : 7(e) 



70 if \9\ < 0j 
1 else 




if 101 < 0, 



else 



(1) 



(2) 
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D 3 : 7(0) = 1 + (70 - 1) exp 
(70 - 1) 



-In 2 



Da : 7 (6») = 1 + 



1 + 



(3) 
(4) 



Velocity profile: spherical 



The symmetrical situation holds for \9\ > 7r/2 (counter-jet): 
7(0)=7(0 + 7r). 

For a given observation angle o b s defined between the 
jet axis and the line of sight, we project on the sky plane 
the surface observed at a given observational time T ^ s . The 
observed flux on the sky plane is related to the intrinsic 
flux in the source rest frame S^.int by the Doppler factor: 
SiAobs — Sv, intf> 3+n - Due to the velocity distribution 7(6), 
each point shell have an intrinsic velocity different in norm 
and direction, and then a different apparent speed as mea- 
sured by the observer. 

The appearance of a moving feature must also be cor- 
rected by time delay effects to the observer. Let us define 
two coordinate frames : IZi (Ozx on Fig. [T]), for which one 
of the axis is aligned with the jet axis, and IZ2 (Oz'x' on 
Fig- HJ, aligned with the observer line of sight, and obtained 
by rotating IZi by an angle 8 bs- If one defines a point M 
with the polar coordinates (r, 8) in IZi , hence the polar coor- 
dinate of M in IZ2 are (r, a = 9 — f? b s )- The observation time 
of M, i.e. the instant when the light emitted at M (when 
the shell just reached it) will be seen by the observer is: 



Tobs(M) 



p(a)c 



P(a)c 



(1- /3(a) cos a) (5) 



where /3(a) is deduced from the velocity distribution 7(6*) 
expressed in 7t2 : 



/3(a) 



7(0 = a + e ohs y 



(6) 



Hence, two points of the jet Mi and Mi will be seen by 
the observer at the same instant if the observation times 
reach the condition: T a b B (Mi) = T b B (M2). Let's choose as 
a reference point the intersection between the propagating 
shell and the jet axis. At a given instant t, this point is at 
distance ro(t) from the origin, and is characterized by the 
Lorentz factor 70. Using Eq. ([5]), the parametric equation of 
the jet surface seen at a given observational time expressed 
in IZ2 is then: 



-(a) = ro(t) 



l(°0 



1 - /3 COS 6> ob s 



1 -£(«)< 



(7) 



No characteristic scale is involved in this equation which 
is self-similar. Thus, the observed surface at a different 
observational time is a simple homothetic transformation of 
the previous one, through the time dependent evolution of 
ro. In the non-relativistic case (/3 <C 1), one can easily check 
that this expression simply yields the intrinsic shape of the 
surface. Figure [2] shows the deformation of the observed 
surface for a purely spherical velocity profile and the 3 
different velocity profiles Di, D2 and D3, using 70 = 3. 
Even for this mildly relativistic regime, the deformation due 
to light propagation delay is strong. Hence, this effect must 
be taken into account if one wants to compute properly 
the observed intensity, projected on the sky plane. For 
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Velocity profile: D2 
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Figure 2. Surfaces seen by an observer at a given arbitrary ob- 
servation time Toba (Eq. [7] with ro = 1), for the non-relativistic 
case (70 = 1.00001, in black), and for the relativistic case (70 = 3, 
in red), for a spherical expanding surface (Upper-left case), and 
the different velocity distributions D\— 3 . The shaded region rep- 
resents the geometrical collimation of the jet (Qj = 20°), and the 
dashed line is the observer line of sight (8 f, s = 30°). 



comparison, the shape of the shell in a non-relativistic case 
(70 — 1 = 10~ 5 ) is shown. In this case, as we remarked, 
the shape is the real one since the light travel effects are 
negligible. 

Knowing the shape of the observed surface, it is now 
possible to project it on the sky plane and to compute the 
intensity profile at each point of the sky. A point M of this 
surface is characterized by its polar coordinate a in IZ2, and 
its radius r(a) given by Eq. (0. The projection of M on the 
sky plane is simply given by: 



x (a) = r(a) sin a? 



(8) 



The apparent velocity (in unit of c) of M is computed 
using the formula for a relativistic point source: 



/3a PP (a) = 



/3(a) sin a x'(a) 



l-/3(a)< 



/3o 



ro V 1 - /3o cos #ob 



(9) 



It is interesting to note that, in our formalism, the apparent 
speed can be expressed in function of the position on the 
sky plane due to the self-similarity of the observed surface. 
The Doppler factor of M is: 



5(a) 



Vi-£(c*) 2 

1 — /3(a) cos a 



(10) 



Assuming a flat spectrum source (n = 0) and a uniform 
intrinsic emissivity, the observed intensity is simply 1(a) = 
Io5(a) 3 . Then the intensity profile is parametrized by these 
two equations: 
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app 



Figure 3. Comparison of the jet apparent velocity for the conical and gaussian profile D\ and D3 and for an observation angle ( 



: 20°, 



a jet opening angle 6j = 15° and a Lorentz factor on the jet axis 70 = 10 . Left: Plots of the two intensity profiles 1(a) (black: conical 
profile, red: gaussian profile), computed following Eq. mil l. Their maxima are not at the same position compared to the line of sight 
(vertical solid line) and are not obtained for the same value of the Lorentz factor. The gray area corresponds to the geometric jet aperture 
i.e. \a — o b s \ 9y Right: Plot of the same intensity profiles but versus the apparent speed j3 app of each point of the shell surface. The 
apparent speed of the whole pattern is then given by the apparent speed of the brightest point of the shell, indicated by the vertical 
dashed lines. In the present example, this apparent speed is equal to 9.85c and 3.4c for the conical and gaussian case respectively. 



x'(a) — ro sin a 
1(a) = I 



Po 

\/l-/3(«) 2 



1 - Po cos 6> obs 

1 — P(a) cos a 
n 3 



/3(a) cos a 



(11) 



2.2 Apparent velocity of a shell 

The above formulae hold for a single emission location. 
However, in real observations, one deals with a complete 
pattern of emission and an observer will derive an apparent 
velocity of the whole pattern. How this apparent velocity 
is defined is not a simple issue and may depend on how 
the different comp o nents are identified and followed. 
iGopal-Krishna et alj (|2007a! ) have adopted a prescrip- 
tion by computing an effective apparent velocity as a 
Doppler-boosted intensity weighted average of the apparent 
vel ocity of each point of the e mitting pattern (cf Eq. (5) 
of IGopal-Krishna et al.1 12007 j ). We think however that 
this prescription does not really reproduce the way by 
which apparent velocities are derived. Usually, observers 
will fit the com ponents by a bell -shape (often gaussian) 
try function (e.g. iLister et al. 200§|) and define the velocity 
by the displacement of the maximum of the fit. Thus we 
have chosen another prescription: we define the effective 
apparent velocity of a component as the velocity of the 
maximum of specific intensity, much like one defines the 
group velocity of a wave packet. 

The difference can be illustrated in a simple case, where 
a conical shell is expanding at a constant velocity. If the 
observer line of sight lies within the cone of expansion, but 
not straight on-axis, he will see an asymmetrical pattern, 
but whose maximum will always be at zero angle with 
the line of sight, since the Doppler factor will always be 
maximal in this direction. An average velocity will thus 



be non-zero, whereas the maximum will not move. To 
illustrate the difference with GK07's prescription, we have 
simulated a conical jet with a constant Lorentz factors 
and a finite opening angle 9j = 15deg. The Lorentz factor 
has been fixed to 76 = 2 or 75 = 5. We note that GK07 
didn't take into account the apparent deformation of the 
surface caused by light travel effects, and have computed 
the averaged velocity for a given proper time and not a 
constant observer time. We thus computed the apparent 
velocity in four different ways : taking the initial GK07's 
prescription without light-travel time correction, a modified 
GK07's prescription taking into account the light travel 
time, our prescription (velocity of the brightest point), and 
the velocity deduced from a gaussian fit of the intensity 
map (with light travel time corrections). The results are 
displayed on Fig. [4] We see that our prescription with 
the brightest point is in complete agreement with what 
would be inferred from a gaussian fit of the intensity map, 
then supporting our own prescription. Inclusion of light 
travel effects in the GK07 prescription produces significant 
differences with their original calculations, especially for 
large Lorentz factor (right panel of Fig. [4]). Interestingly, 
this "modified" prescription looks rather similar to our 
results. 

Due to the axial symmetry hypothesis on the jet geome- 
try, the problem of determining the position of the brightest 
point of a tridimensional surface projected on the sky plane 
can be treated in a bi-dimensional approach. Indeed, the 
maximum of intensity is necessarily in the plane defined by 
the jet axis and the observer line of sight, i.e. the plane of 
Fig. [1] The position of the maximum is computed numer- 
ically , using a "Golden Section search" algorithm (|Kiefeil 
Il953l ). As an example, Fig. [3] shows a "visual" comparison 
of the different steps of the computation of the jet apparent 
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Figure 4. Comparison of four different prescriptions for the apparent velocity of a thin shell. The left panel is computed with 70 = 2, 
and the right panel with 70 = 5. Solid line: this paper's prescription. Dot-dashed line: original GK07's prescription with an averaged 
velocity, without light-travel effects; Dashed line: modified prescription of GK07 including light travel effects. Red dots: results of the 
gaussian fit with a simulated image. 
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Figure 5. Apparent velocity in unit of c (black open circle), and 
Dopplcr factor (black open diamond), as a function of # b s , com- 
puted with our conical jet model (velocity profile D\) assuming 
Oj = 0.001 and 70 = 10. We also display the theoretical curves 
for the apparent velocity (red line) and Dopplcr factor (red dot- 
ted line) of a point source of same Lorentz factor 70 = 10. The 
agreement is excellent. 

speed for the conical and gaussian profile D\ and D3. The 
two intensity profiles 1(a), computed following Eq. ((TTJ , are 
plotted in the left panel of Fig. [3] Their maxima are not at 
the same position compared to the line of sight and are not 
obtained for the same value of the Lorentz factor. This re- 
sults in a different apparent velocity of the whole pattern on 
the sky plane as shown in the right panel of Fig. [3] In con- 
sequence the velocities of the maximum of intensity, which 
correspond to the apparent velocities of the jet in our for- 
malism, are different of the order of 9.85c for the D\ profile, 
and 3.4c for the D2 profile. 

2.3 Test case: perfectly collimated jet 

We have tested our method to compute the appearance of 
a relativistic jets (apparent velocity and Doppler factor) in 
the case of a very narrow collimated jet. The geometrical 
angle being much smaller than the radiation cone, the 
result should be very close to the point source case. The 
agreement is indeed excellent as shown in Fig. [5] our model 
of very well collimated jet mimicking perfectly the behavior 



of a point like source. 



3 RESULTS 

We are now able to compute the apparent velocity of a jet 
with a given velocity profile [Oj, 70, Di), and for a given 
observation angle o ^ B . By iterating the process for differ- 
ent o bs, we are able to compute the relation between the 
observation angle and the apparent speed of the whole pat- 
tern, and to determine the maximum apparent velocity that 
can be reached for a given jet configuration, as well as the 
associated Doppler factor, for all possible orientations. 

3.1 Effect of the jet velocity profile 

Figure [5] shows the evolution of the jet apparent velocity and 
the Doppler factor as a function of the observation angle 8 h s 
for the four velocity profiles Di_4 defined in Eqs.[TJ-|4l We 
assume a geometrical collimation 9j — 15° and a Lorentz 
factor on the jet axis fixed to 70 = 10. As a comparison, we 
also display the theoretical curves of the apparent velocity 
(red line) and Doppler factor (red dashed line) expected in 
the case of a perfectly collimated jet with the same Lorentz 
factor. 

3.1.1 Conical jet 

In the case of the conical velocity profile Di (upper-left 
corner of Fig. , and as we noted previously, the apparent 
velocity is equal to for o h s ^ 0j, while the Doppler factor 
is maximal (S = 270). For a higher observation angle, the 
apparent velocity and the Doppler factor follow the same 
behavior as the theoretical expressions, but shifted by an 
angle 0j. 

These results can be understood rather easily: while the 
observer is looking inside the jet cone (0 o bs ^ Oj), the ob- 
served image is dominated by the part of the shell that points 
directly toward the observer. Then, the measured apparent 
velocity is null, and the Doppler factor is maximal (270). 
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Velocity profile D 1 Velocity profile D 2 




Figure 6. Jet apparent velocity (black empty circles) in unit of c and Doppler factor (black empty diamonds), as a function of the 
observation angle # b s , computed by our model for different jet velocity profiles. Upper- left corner: velocity distribution D\ (Eq. [TJ, 
upper-right corner: velocity distribution D2 (Eq. [2J, lower-left corner: velocity distribution D3 (Eq. [3}, lower-right corner: velocity 
distribution D4 (Eq. [4j. The jet opening angle is 8j = 15 °, and the Lorentz factor on the jet axis is 70 = 10. For comparison, we have 
also plotted the theoretical curves for the apparent velocity (red line) and the Doppler factor (red dashed line) of a point source with a 
Lorentz factor of 10. 



In this case, the image seen by the observer will be a dis- 
torted circle, expanding with time homothetically around 
the brightest point of the image. 

As soon as the line of sight is out of the cone (8 bs > 0j), 
the image is dominated by the emission coming from the 
edge of the conical jet. Then, the situation becomes similar 
to a point-like source, seen under the angle 9 b s — Oj- 
The evolution of the apparent velocity /3 ap p(#obs) and 
the Doppler factor 5(9 bs) is then the same as for the 
theoretical one in case of a perfectly collimated jet, shifted 
by 8j , the external ridge playing the role of a thin narrow jet. 

For this kind of velocity profile, it is always possible to 
observe high super-luminal motion if the jet is seen close to 
the edge of the cone. For each cone angle 8j, the maximal 
apparent speed reachable is /3 app , max ~ 70 for an observation 
angle r? b s ^ 9j + To" 1 - The Doppler factor is then S ~ 70. 
But more interestingly, a wide range of observation angles 
can be associated to an absence of apparent motion, while 
the Doppler factor is maximum. This situation is unlikely if 
one consider a perfectly collimated jet, but become highly 
probable if the jet opening angle is large. 



3.1.2 Jet with a continuous velocity profile 

Profile D2 is obtained from D\ by adding a power-law 
decrease of the Lorentz factor outside the jet cone (Eq. [2]), 



so that the jet velocity smoothly decreases to zero at 
large 9. The corresponding curves of the shell apparent 
velocity and Doppler factor as a function of the observation 
angle (upper-right corner of Fig. [6]) are the result of a 
complicated convolution of the Lorentz factor profile, with 
the theoretical formula for a point-source. Similar to the 
Di profile, the apparent velocity is null when the observer 
looks inside the cone (0 o bs < 9j). As soon as the observation 
angle is larger than the collimation angle (# b s > 9j), 
the brightest point of the jet moves out from the line of 
sight, yielding to a non-zero apparent velocity. However, 
for a given observation angle 9 a i, B , the apparent velocity is 
lower than in the case of D\. This is due to 2 reasons: 1) 
the lower intrinsic Lorentz factor 7 at the position of the 
brightest point on the shell, compared to the Lorentz factor 
value 70 on the jet axis, and 2) this brightest point is very 
close to the line of sight with a < 7 _1 . This last point is 
confirmed by the high value of the Doppler factor (S > /3 app ). 

The lower part of Fig. [S] shows the results of the model 
for a gaussian velocity profile (-D3, Eq. [3} and a lorentzian 
profile (D4, Eq. [4j . For both distributions the apparent ve- 
locity is never null for any observation angle but # b s = , 
the brightest point of the shell being always slightly shifted 
out the line of sight. However, we emphasize that the appar- 
ent velocity is significantly smaller compared to the point- 
like source, for the same two reasons explained before. As 
for the D2 velocity profile, this is confirmed by the high 
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values of Doppler factor: S > /3 a 



v e a 



We can observe 



in Fig. [6] that the maximal apparent speed is reached for an 
observation angle close to the jet opening angle 9 b s ~ 9j±e. 



3.1.3 Analytical approximation for the apparent velocity 

The previous analysis seems to indicate that the maximal 
apparent speed that can be observed with a given velocity 
profile is linked to the angular velocity gradient inside the 
jet. Indeed, /3 app is maximal with D\, for which the velocity 
discontinuity at the edge of the cone induces a large gradi- 
ent. The lowest /3 apP are obtained with the Lorentzian pro- 
file (-D4) which is the smoothest profile (for angles greater 
than 8j) among the four. Furthermore, /3 app is null when 
the observer points towards a constant velocity region in 
the jet (see the case of D\ and D2 with 8 ^ s < 8j, upper 
part of Fig. [6]). The exact analytical study of the expres- 
sion which gives the apparent velocity /3 app (6'obs) is pretty 
involved. However, a Taylor development at the first order 
of the Doppler factor as a function of the angle a in IZ2 can 
explain this general behavior. 

Indeed, for high Lorentz factors (7 3> 1) and small ob- 
servation angles (a < j~ ), the Doppler factor can be writ- 
ten as: 

2 7 



S(a,-y) 



(12) 



1 + 7 2 Qf 2 

In our model, we make the assumption that the intrinsic jet 
emissivity is constant. Hence, the position of maximum ob- 
served intensity corresponds to the position of the maximum 
Doppler factor. Thus the logarithmic derivative of the previ- 
ous equation vanishes at the position of maximum observed 
intensity: 



d<5 
T 



d7 27Q 2 d7 + 2cry 2 da 



= 



(13) 



7 1 + 7 2 a 2 

After some calculation, this equation can be reduced to a 
simple quadratic equation for a: 



9 7 1 , d7 

a + 2-a = with 7 = -r- 

7 7 2 da 

The two solutions of this equation are: 



a± 



1 ± 



1 + ^ 

7 4 



(14) 



(15) 



We will consider two cases, following the sharpness of 
the velocity profile. A smooth profile will be defined by \ j\ <C 
7 2 , meaning that the decrease of 7 takes place over an angle 
7 1 

interval Aa ~ — 3> — . 

M 7 

the previous expression can be linearized and we obtain the 
solutions: 



In this case, the square root in 



a+ 



+ 



27 s 



and 



7 
2 7 3 



(16) 



We can easily check that the first solution corresponds to a 
minimum of Doppler factor because |a+| 3> I/7. Then, the 
maximum of intensity is reached for the second solution a_. 
The apparent velocity computed at this position is : 

/3sin|a_| 27 2 |a_| 



(17) 



1— /?cosa_ l + 7 2 a 2 _ 
Now if |a_ I < I/7, and using Eq. (|16|l . this reduces to 



/3 app ~ ill ~ Aa" 
7 



(18) 



We can also check that an absence of velocity gradient 
in the jet (7 = 0) implies that the maximum of intensity 
is aligned with the line of sight (a_ = 0), resulting in the 
absence of apparent velocity. This is a very interesting 
conclusion: for a smooth angular variation of the Lorentz 
factor, the apparent velocity of the jet is linked to the width 
of the angular velocity profile, and not on the absolute 
value of the maximal Lorentz factor. 

In the case of a sharp variation, defined by I7I 3> 7 2 , 
the linearization of the square root yields a different result. 
We obtain two symmetrical solutions, independent of the 
velocity gradient: 



a± 



±1 

7 



(19) 



This applies to the case of the conical jet profile D 1 , for 
which there is a discontinuity in the velocity distribution. 
This discontinuity implies a infinite gradient at the edge 
of the jet. We can check on Fig. [5] that the maximal jet 
apparent velocity (/? app = 70) is reached for 0=1/7 i- e - f° r 
an observation angle 8 h s = 6j 1 + 1/7. In this case, the range 
of apparent velocities is the same as for a point source (see 
Fig. [SJ) as discussed in Sect. 13.1.11 



3.2 Maximal apparent velocity 

3. 2. 1 Effect of the jet opening angle 

For a given jet configuration {Di, 70, 8j), we now define 
/3 app , ma x, the maximal apparent jet velocity that can be 
obtained by varying the observation angle # bs- In order 
to study the influence of the geometrical collimation, we 
have computed the evolution of /3 app , ma x for the different 
velocity profile D1-4, as a function of the jet opening angle 
8j. The results are presented in Fig. [7J where we have 
plotted /3 apP , ma x(#j) and the associated Doppler factor 
5p a pp,max(9j) f° r eacn velocity profile and assuming a 
Lorentz factor 70 = 10 on the jet axis. 

For the profile D± (conical profile), the maximal appar- 
ent velocity is constant, /3 app , max ~ 70. The corresponding 

ax . As 



Doppler factor 



is then of the order of /3 a 



already mentioned, this profile can be seen as a point source 
when the observer looks outside the cone (# bs > 8j). Then, 
it is always possible to find Sobs that matches the conditions 
to obtain the maximal theoretical value /3 app = 70- 

The three other profiles D2-4 have all a similar behavior 
which is however very different from the Di case. For small 
opening angles (8j <C "f^ 1 ), /3 app . max is close to the theoret- 
ical maximum because in that case the jet is very well colli- 
mated and can be seen as a point-like source. When 8j starts 
to increase, the value of /3 app , max decreases rapidly. Depend- 
ing on the velocity profile, the corresponding Doppler factor 
increases or decreases, but remains much higher than the 
value of /3 app , max . Hence, taking into account the effect of 
the jet opening angle allows to obtained small apparent ve- 
locities but with higher Doppler factors than in the homoki- 
netic case. 



8 T. Boutelier et al. 



Maximal apparent velocity 






Yo=10 






Profile 


Di : 


—A— 


Profile 


D 2 : 


• 4- 


Profile 


d 3 : 




Profile 


d 4 : 



.aims 



15 
^10 



Doppler factor 

' I ' ' ' H 

4 A- A-A-A A A- A-A ' 



T 



OOOOK X- X -X X-X- X-X- X- X -X" 



20 40 

9: 



60 



5 r 

oEl 





20 40 

9: 



60 



30 
25 

S20 

E 
£5 

^10 



Maximal apparent velocity 



4 





Yo = 30 






Profile 


Di : 


—A— 


Profile 


D 2 : 


• +• 


Profile 


D 3 : 


■ X 


Profile 


D 4 : 



OEL-i ■ ■ I ■ ■ ■ 1 ■ ■ ■ +- 



60 F 
50 
40 
=^30 
20 
10 



Doppler factor 

a-aI A A-a-^a-a' A £ 



A 



^ ^30 ' 4, 

- on - 



DOooooe<x»-x x-x * x-x X-X X- X- 



" + ++ + . 



OEL 



I i i i I i i i £ 



20 



e 



40 



60 



20 



e 



40 



60 



Figure 7. Left: Evolution of the maximal apparent velocity (/3 a pp,max) as a function of the jet opening angle for the different velocity 
profiles D1—4. Blue circles: profile D\, violet triangles: profile D2, green crosses: profile D3, salmon crosses: profile D4. For all profiles, 
we have set 70 = fO (top) or 70 = 30 (bottom). Right: Same as left, but for the Doppler factor. 



3.2.2 Empirical estimate 0} the jet collimation 

We have computed, for each velocity profile, surfaces of the 
maximum apparent velocity /3 apP , max and the associated 
Doppler factor 5p max in the (9j,yo) plane. Plots reported 
in Fig. [7] then correspond to a section of these surfaces at a 
constant 70. 

These surfaces are plotted in Fig. [S] with contours 



of /3 a 



and 5b 



in solid lines. In the case of 



the conical profile (top panel of Fig. (8|, the maximum 
apparent velocity follows the theoretical expectations of 
a relativistic point source (see previous section) with 
/3ap P ,max(#j,7o) — 7o and a corresponding Doppler factor 
Sp apPimax — 70, both parameters being independent of 
the jet opening angle 8j. The other velocity profiles give 
completely different dependencies on 9j but also on 70 
(see the 3 lower panels of Fig. [5]), the maximum apparent 
velocity becoming strongly dependent on the jet opening 
angle, while the corresponding Doppler factor is very 
sensitive to the velocity profile. It confirms the affirmations 
that we make in the previous sections, i.e. that the jet 
opening angle decreases dramatically the apparent velocity, 
even for high jet Lorentz factor. And the larger the jet 
opening angle, the smaller the value of /3 a pp,max compared 
to the point-source case. 

But more interestingly, Fig. [8] allows to determine an 
empiric relation between the jet apparent velocity an its 



collimation angle for the continuous velocity profiles D2-4: 
for instance, the collimation must be better than 9j < 5 , 
in order to observe an apparent velocity of at least 10c (i.e. 
fiapp,max > 10). This limit on the collimation can be ex- 
pressed as 5 ~ 0.1 rad ~ Pap P) max an d can then be re- 
written as: 



9j < PI 



(20) 



This relation remains true by a factor of a few for the pro- 
files D2-4 but appears to better work for continuous ve- 
locity profiles like D3 and D4. The abrupt variation of 7 in 
Di helps in reaching large apparent velocity for jet opening 
angle a bit larger than Pa P \, max ■ However, since we expect 
continuous profiles to be more realistic, Eq. (|20[> should be 
relatively general. Hence, the measurement of the apparent 
velocity of a relativistic jet can provide an upper limit on 
its collimation angle. 



3.3 Effect of an angular dependent emissivity 

In all the previous sections, we assumed a constant specific 
intensity in the shell frame, in all directions. The variation 
of apparent brightness is thus entirely due to varying 
Doppler factors. This is obviously a very crude assumption 
since the physical reasons causing a variation of jet velocity 
are also likely to produce variations of non thermal acceler- 
ation processes and thus of non thermal radiation emission. 



The influence of the jet opening angle on the appearance of relativistic jets 9 



Profile r>i 

Maximal apparent velocity Corresponding Doppler factor 




Profile D 2 

Maximal apparent velocity Corresponding Doppler factor 




Profile r> 3 

Maximal apparent velocity Corresponding Doppler factor 




Profile D A 

Maximal apparent velocity Corresponding Doppler factor 




Figure 8. Contour plots of the maximal apparent velocity (/3 a pp,max on the left) and the associated Doppler factor (on the right), as 
a function of the geometrical collimation of the jet 6j and of the Lorentz factor on the jet axis 70, for the 4 different angular velocity 
profiles (from top to bottom) D\, D2, D3 and D4 (Eg. 1114) . 
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Figure 9. Evolution of the apparent velocity as a function of the 
observation angle for a constant intensity profile (thick solid line) 
or a angular dependent (gaussian in the present case) intensity 
profile (thick dashed line). We use a gaussian profile (-D3) for the 
Lorenz factor. For comparison the case of a conical profile (Di) 
and constant intensity is plotted in thin line. The jet opening 
angle 6j = 15 and 70 = 5 



Considering a possible variation of local emissivity with 
the angle will add another complexity to the problem. We 
can remark however that for a given direction of the line 
of sight, a variable specific intensity can be considered as 
equivalent to a modification of the Doppler factor, since 
it is always possible to find another velocity producing 
the same apparent specific intensity at each position on 
the sky. Letting the specific intensity free to vary won't 
probably change much the qualitative conclusions about 
the limitation of apparent velocities and Doppler factors. 

We illustrate the influence of a variable intensity with 
the D3 distribution (gaussian profile) of the Lorenz factor 
and two different assumptions on the specific intensity. The 
results are displayed in Fig. [9] for 70 = 5. We chose either 
a constant profile of the specific intensity (thick solid line) 
as computed previously, or a gaussian profile (thick dashed 
line). As already noticed, a constant intensity and Lorentz 
factor within a finite angle 8j (conical profile) will always 
give a vanishing apparent velocity for o b s < Oj since the 
maximal Doppler factor will always been along the line of 
sight, giving a null projection on the sky. For o t s > Oj the 
image is dominated by the emission coming from the edge 
of the conical jet and the situation becomes similar to a 
point-like source. This is exemplified by the thin solid line 
in Fig. [9] Things are different with a variable Lorenz factor 
and variable intensity, since the maximum is no more strictly 
on the light of sight (it is determined by some compromise 
between the variation of the Doppler factor and the specific 
intensity). The apparent velocity will also be not zero, since 
it is directly related to the angular separation between the 
maximum and the line of sight. However, we note that the 
resulting pattern is not strongly modified by the inclusion of 
an angular dependency of the specific intensity profile and 
most of the effects are already catched by using an angle- 
dependent profile for the Lorenz factor. Since we don't have 
anyway a good knowledge of the real distribution of the 
Lorentz factors, the inclusion of a angular dependency of 
the emissivity wouldn't really add any new effects. 



4 DISCUSSION 

All these simulations show that the apparent velocity of a 
jet is highly dependent on the geometrical collimation and 
the velocity profile inside the jet. The apparent velocity 
decreases dramatically compared to the homokinetic case 
as soon as the jet collimation angle is larger than a few 
degrees. These resu lt s are in a g reemen t with those of 
iGopal-Krishna et all (|2004 120061 , l2007al lbh. although the 
inclusion of some effects neglected by these authors modify 
the precise distribution. Our choice of another prescription 
to define the apparent velocity, as the velocity of the 
maximum of the specific intensity of the observed pattern, 
gives rise also to different values: for instance, as we showed 
in Sect. 12.21 the number of objects having a null apparent 
velocity is finite in our case for a conical, single velocity jet, 
whereas it is true only exactly on-axis with the prescription 
of GK07 where the apparent superluminal velocity is 
obtained by averaging the different apparent velocities of 
the observed pattern. 

As shown in Fig. [6] the distribution of apparent veloci- 
ties is significantly different from that deduced from a single 
velocity Doppler boosted jet. This can strongly modify the 
statistical distribution of apparent velocities and thus the 
bulk Lorentz factors deduced from these statistical studies. 

In the case of a smooth decrease of the bulk Lorentz 
factor with the angle to the LOS, we have shown that the 
apparent velocities are much more indicative of the geo- 
metrical shapes of the jet than of the real bulk Lorentz 
factors. Thus the bulk Lorentz factors deduced from ra- 
dio observations and statistical studies can be very dif- 
ferent from those deduced from radiative constraints for 
high energy emission, and will give generally much lower 
values. This can obviously help solving the discrepancies 
between the two estimates, since Lorentz factors deduced 
from gamma-rays transparency arguments routinely exceeds 
several tens for TeV blazars, alt hough superluminal mo- 



tion are seldom 



Piner &i Edwards! 



observ ed 
12004 



( e.g. lEdwards fc Pinerl |2002|: 



Tavecchio. Maraschi. &: Ghiselliiul 



19981 ; ISauge fc Henrill2004 ). Note however that the discrep- 
ancy can be partly reduced w ith inhomogeneous models of 
jets (e.g. lHenri fc Sauge|[20of3 ). 

A sharp decrease of the Lorentz factor outside a cone of 
almost constant Lorentz factors will give a different appear- 
ance: in this case, the apparent superluminal motion will 
be very small, maybe undetectable, inside the collimation 
cone, although the Doppler factor is maximal. But a small 
fraction of objects, seen at a small angle from the edge 
of the cone, can exhibit large superluminal motions. Both 
zero apparent velocities and highest apparent velocities will 
be overrepresented in statistical samples compared to the 
homokinetic hypothesis. These possibilities should of course 
be carefully taken into account in statistical studies. 

Contrary to apparent superluminal velocities, Fig. [8] 
shows that the effective Doppler factor strongly depends 
on the precise velocity profile inside the jet. This can be 
understood by remarking that, in the case of a relativistic 
point source (Fig. 0, when the apparent velocity is max- 
imal, near 6 t B — Jq , the Doppler factor is varying very 
rapidly. Thus, small differences in the precise angle of the 
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emitting regions contributing the most to the intensity can 
lead to significative differences in the Doppler factor, even 
if the apparent velocities are comparable. For instance, in- 
creasing the opening angle of the jet will generally decrease 
the apparent velocity, but can decrease or increase the effec- 
tive Doppler factor (see Fig. 0. 

In consequence, in the case of open jets, it is difficult to 
derive a statistics of Doppler factors from the observation 
of apparent superluminal motions without knowing the 
precise angular distribution of Lorentz factors. 

These considerations may give a hint to solve several 
discrepancies that arise when comparing superluminal mo- 
tions and other constraints derived from radiative models 
or statistical arguments. Namely many TeV blazars show 
little or no apparent motion, whereas radiative constraints 
seem to be compatible only with very high bulk Lorentz 
factors. As most of TeV blazars are BL Lacs, which are 
thought to be beamed counterparts of weakly collimated 
FRI galaxies, it is tempting to explain this discrepancy 
by a rather large opening angle. It is interesting to note 
that very open jets are commonly invoked in the context of 
gamma-ray bursts, where the usually assumed bulk Lorentz 
factors are much higher (a few hundreds) than 6J 1 (se e 
however the "cannonball" model, iDar & De Ruiulall2000h . 
The visual appearance of a gamma-ray burst in the earlier 
phase would be thus dominated by the part of the jet 
traveling at zero angle with respect to the line of sight, i.e. 
with no superluminal motion at all. 

It should be stressed that we studied only a variation 
of 7 with the angle of ejection, but that other possibilities 
exist. One could also imagine that several (and may be a 
whole distribution) of Lorentz factors exist along a partic- 
ular direction, even for a perfectly collimated jet. This is 
indee d necessary for all "shock in jets" models (see IZensusI 
Il997l for a review), where layers at different velocities are 
supposed to be ejected and collide through inter nal shocks. 
So called "spine in jets" and "two-flows" models (|Ghisellinil 
120051 : iHenri fe Pelletierill99ll ) also assume different Lorentz 
factors along a single direction. We didn't consider either the 
possible variation of the bulk Lorentz factors with the dis- 
tance. All these factors would of course complicate the pic- 
ture, since they would produce different apparent velocities 
and bulk Lorentz factors along the jet, and most probably 
varying following the observed wavelengths. 



5 CONCLUSION 

We have studied the influence of the geometrical opening of 
relativistic jets on their apparent velocities and associated 
Doppler factors. We parametrized this opening by various 
distributions of Lorentz factor as a function of the angle of 
observation. We have investigated the appearance of a thin 
shell ejected and traveling along the flow. We propose a new 
criterion to define the apparent velocity of the jet pattern 
as the velocity of the point having the highest observed 
specific intensity. The effective Doppler factor is then the 
Doppler factor associated with this point. 

Studies of different configurations shows that the 



variation of the maximal apparent velocities as a function 
of the viewing angle can be strongly modified, especially 
when the Lorentz factor angular distribution is decreasing 
smoothly. In this case, the maximal apparent velocity is 
essentially determined by the typical angular interval on 
which the decrease occurs (which is also the opening angle 
of the jet for reasonable configurations), and not by the 
bulk Lorentz factor. The associated Doppler factor is a 
sensitive function of the precise shape of the distribution. 
This could help understanding the discrepancies between 
various estimates of bulk Lorentz factors in extragalactic 
jets. 

These findings could have significant impacts on sta- 
tistical studies trying to estimate Doppler factors from 
statistical comparisons between beamed and unbeamed 
objec ts, in the frame of s o-called "unification" models 
(e.g. lUrrv fe Padovanil Il995l ). The opening of jets modi- 
fies strongly the angular distribution of observed effective 
Doppler factors, and thus the statistical distribution for ran- 
domly oriented jets. It could thus significantly alter the con- 
clusions of such studies. Further work will be done to inves- 
tigate this issue. 
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